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The spectrum of a system with multiple channels composed of two hadrons with nonzero total 
momentum is determined in a finite volume with periodic boundary conditions using effective field 
theory methods. The results presented are accurate up to exponentially suppressed corrections in 
the volume. The formalism allows one to determine the phase shifts and mixing parameters of 
Tvn — KK isosinglet coupled channels directly from Lattice Quantum Chromodynamics. We show 
that the extension to more than two channels is straightforward and present the result for the three 
channels. From the energy quantization condition, the volume dependence of electroweak matrix 
elements of two-hadron processes is extracted, for both relativistic and non-relativistic systems. In 
the NR case, we pay close attention to processes that mix the ^So — ^Si two-nucleon states, e.g. 
proton-proton fusion {pp d + e'^ + u^) , and show how to determine the transition amplitude of 
such processes directly from lattice QCD. 

Scattering processes in hadronic physics provide useful information about the properties of particles and their 
interactions. Some of these processes are well investigated in experiments with reliable precision. However, there are 
interesting two-body hadronic processes whose experimental determinations continue to pose challenges. They mainly 
include two-body hadronic scatterings near or above the kinematic threshold with the possibility of the occurrence 
of resonances. Here we discuss two pertinent cases in Quantum Chromodynamics (QCD), the first of which is the 
scalar sector, whose nature is still puzzling (See for example [1 and references therein). While some phenomenological 
models suggest the scalar resonances to be tetraquark states (as first proposed by Jaffe [2j), others propose these to 
be weakly bound mesonic molecular states. The most famous of which are the flavorless ao(980) and /o(980), which 
are considered to be candidates for a KK molecular states [3] , [3] , [5] . In order to shine a light on the nature of these 
isosinglet resonances, it would be necessary to perform model-independent multichannel calculations including the 
{tttt, tttttttt, KK, 7777} scattering states directly from the underlying theory of QCD. In the baryonic sector, observations 
of the strong attractive nature of the isosinglet KN scattering channel led to the postulation of kaon condensation 
in dense nuclear matter (6 . However, extracting KN scattering parameters is a rather challenging task, due to the 
presence of the Ett scattering channel and the A (1405) resonance bellow the KN threshold, see for example Ref. [7]. 
Previous chiral perturbation theory (^PT) calculations have found inconsistency between experimental determination 
of the KN scattering length from scattering data and kaonic hydrogen level shifts ^8], [9], [10], [11], [12], but as with 
any low energy effective field theory (EFT) calculation, there are unaccountable systematic errors associated with 
the large number of unknown low energy coefficients (LECs) needed to perform accurate calculations of multichannel 
processes. 

In addition to these strongly coupled scattering processes, there are weak processes involving multi-hadron states 
that require further investigation. For instance. Lattice QCD (LQCD) calculations have recently shown further 
evidence of S^S^ and AA (_ff-dibaryon) shallow bound states [13], [14], [15]. This will certainly reignite ex- 
perimental searches for evidences of these states. Among the possible weak decays of the iJ-dibaryon include 
H — >■ {nA,nYj'^,pYj~,nn) \W . In hyper-nuclear physics, there has been much interest in a definitive determination 
of the contribution of non-mesonic weak decays (AiV — > NN, ANN — NNN) to the overall decay of hyper-nuclei. 
In particular, as discussed in |17| (and references therein), there has been a long standing puzzle regarding the theo- 
retical underestimation of the ratio of the decay widths r(An — > nn)/T{Ap — > np) as compared to the experimental 
value. Certainly a great deal of progress has been made by meson exchange models in order to close this gap, but 
a model-independent calculation directly form QCD would give further insight into the mechanism of these decays. 
These two cases illustrate processes where it is necessary to evaluate weak matrix elements between multi-hadronic 
states. 

Currently, LQCD provides the most reliable option for performing calculations of low energy QCD observables. 
LQCD calculations are necessarily performed in a Eucledian and finite spacetime volume. Although the former 
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forbids one to calculate the physical scattering amplitudes from their Euclidean counterparts away from the kinematic 
threshold due to the Maiani- Testa theorem [18], the latter is proven to be a useful tool in extracting the physical 
scattering quantities from lattice calculations. In his prominent work, Liischer showed how one can obtain the infinite 
volume scattering phase shifts by calculating energy levels of two-body scattering states in the finite volume (TH] , [50] ■ 
The Liischer method which was later generalized to the moving frames in Ref. pT^E5] , is only applicable to scattering 
processes below the inelastic threshold. Therefore it cannot be used near the inelastic threshold where new channels 
open up, and a generalized formalism has to be developed to deal with the coupled(multi)-channel processes. A direct 
calculation of the near threshold scattering quantities using LQCD can lead to the identification of resonances in 
QCD such as those discussed above, and provide reliable predictions for their masses and their decay width. One 
such generalization was developed by Liu et al. in the context of quantum mechanical two-body scattering |24j . 
[25j . There, the authors have been able to deduce the relation between the infinite volume coupled channel S-matrix 
elements and the energy shifts of the scattering particles in the finite volume by solving the coupled Schrodinger 
equation both in infinite volume and on a torus. The idea is that as long as the exponential volume corrections are 
sufficiently small, the polarization effects, as well as other field theory effects, are negligible. Therefore after replacing 
the non-relativistic dispersion relations with their relativistic counterparts, the quantum mechanical result of Liu 
et al. [21], [IS] is speculated to be applicable to the massive field theory. In another approach, Lage et al. considered 
a two-channel Lippman-Schwinger equation in a non-relativistic effective field theory (NR EFT). They presented the 
mechanism for obtaining the KN scattering length, and studying the nature of the A (1405) resonance from LQCD 
[26] . Later on, Bernard et al. generalized this method to the relativistic EFT which would be applicable for coupled 
KK — TTTT channels [T]. Unitarized chiral perturbation theory (UCHPT) provides another tool to study a variety of 
resonances in the coupled channel scatterings. This method uses the Bethe-Salpeter equation for a coupled- channel 
system to dynamically generate the resonances in both light meson sector and meson-baryon sector in the infinite 
volume, see for example [27], [5S], [21], [SO]. When applied in the finite volume, the volume-dependent discrete energy 
spectrum can be produced, and by fitting the parameters of the chiral potential to the measured energy spectrum on 
the lattice, the resonances can be located by solving the scattering equations in the infinite volume. This method has 
been recently used to study the resonances /o(600), /o(980) and ao(980) in [3T], [21], A(1405) in [Mj, ai(1260) in [31], 
Ac(2595) in [3S], and _Ds*o(2317) in jM] m the finite volume. One should note that in contrast to the single channel 
scattering in the finite volume, the coupled-channel scattering requires determination of three independent scattering 
parameters which would require in the very least three measurements of the energy levels in the finite volume. As 
proposed in [1], |31j . one can impose twisted boundary condition in the lattice calculation to be able to increase the 
number of measurements by varying the twist angle and further constrain the scattering parameters. Another tool 
to circumvent this problem is the use of asymmetric lattices as is investigated in [T] , |31j , [33] . Alternatively, one can 
perform calculations with different boost momenta [32], [34| . 

The goal of this paper is two-fold. First, we present a model-independent fully relativistic framework for determining 
the finite volume (FV) coupled-channel spectrum in a moving frame. Secondly, we show how to extract current- 
operator matrix elements in the two hadron sector directly from LQCD, for both relativistic and non-relativistic 
processes. These two are small stepping-stones towards one of the overarching goals of hadronic physics, which is to 
determine properties of multi-hadron system directly from the underlying fundamental theory of QCD. 

This paper is structured as follows. In the first section we present the result from a scalar field theory model, which 
illustrates all the features of the problem at hand and allows us to derive the quantization condition (QC) using a 
diagrammatic expansion. Although, in deriving the QC for the toy model we make a series of approximations, it is 
shown that this result is in perfect agreement with the exact QC which is obtained from the generalization of the 
work by Kim et al. |22j . |37j . In section IB we derive the general form of the finite volume quantization condition 
for N arbitrarily strongly coupled two-body states. This result has been independently derived and confirmed in a 
parallel work by Hansen and Sharpe [55] • In this paper, most of the emphasis will be placed on the N=2 case, but 
the result for the N=3 case will be also explicitly shown. 

After developing the finite volume coupled-channel formalism, we extend our work to be able to determine elec- 
troweak matrix elements in the two hadron sector. The formalism for extracting the physical transition amplitude for 
K — ^ TTTT from the finite volume matrix elements of the weak Hamiltonian in the finite volume, has been developed by 
Lellouch and Liischer [31]. This former quantity has been shown to be proportional to the latter at the leading order 
in the weak coupling, and the proportionality factor, the so-called LL factor, is shown to be related to the derivative 
of the TTTT scattering phase shift. The generalization of the LL factor to the moving frame is given in [22], [37) . 
Ultimately, it would be desirable to calculate electroweak matrix elements of states containing an arbitrary number 
of hadrons, but for the time being we restrict ourselves to the two body sector. Nevertheless, the two body sector 
encompasses a large number of interesting processes including those discussed above. Additionally, as is discussed in 
[40j . the weak disintegration of deuteron in processes such as vd — > vpn and z/gd — > nne'^ is of great importance in 
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neutrino experiments. These processes entail an electroweak mixing between ^5*0 Si two-nucleon channels, and for 
energies bellow the pion-production threshold the dynamics of this system can be described using an EFT without 
pions {EFT {]()) [41], [42], [43]. One hopes to get the physical weak transition amplitude in these processes from a 
LQCD calculation. By applying two alternative methods, degenerate perturbation theory and the density of states 
method, we give a generalization of the LL factor for 2 — >■ 2 weak processes in a moving frame. This includes the 
derivation of the relativistic LL factor in the mesonic sector as well as its non-relativistic extension for the baryonic 
sector. Finally the effect of one-body weak interactions on the finite volume quantization condition and on the LL 
factor is discussed. 



I. MESON-MESON COUPLED CHANNELS 



The goal of this section is to present the quantization condition for N coupled-channel system in a moving frame. 
We present two independent ways to obtain the quantization condition desired. In section |I A[ we present a toy 
model that illustrates the features of the problem at hand and allows us to determine the scattering matrix using a 
diagrammatic expansion. In order to derive the quantization using this approach, it is convenient to consider the case 
where the parameter responsible for mixing the two channels is small. However, as shown in section |I B[ the result 
presented is exact for an arbitrary mixing parameter. 

In section |I B[ we show how to obtain the quantization condition for N arbitrarily strongly coupled two-body 
channels in a moving frame by solving the Bethe-Salpeter equation. This case is the most relevant if we want to 
consider systems coupled via QCD interactions, e.g. the / = system nn (KK, rjrj) — >■ tttt F] 



A. Toy problem: two weakly coupled two-meson channels in the moving frame 

In this section, we consider a two-meson coupled system with J=:0 angular momentum. The two channels will be 
labeled / and //. In general the four mesons can have different masses and quantum numbers, but for the time 
being we restrict ourselves to the case where the two mesons in each channel are identical. Using the "barred" 
parameterization [44) . the S-matrix describing this system can be written as 

f e*2'^^cos2e ie*(''^+*") sin 2e\ 

where 6i and Su are the phase shifts corresponding to the scattering in channels / and // respectively, and e is a 
parameter which characterizes the mixing between the channels. The subscript 2 on 5 denotes the number of coupled 
channels. 

At energies below the four-meson threshold, the dynamics of such system can be described by a simple scalar 
effective field theory (EFT) 

^ = - + [-2 ^ j [g/2 cn) + 

where (f>i is the meson annihilation operator for the i*'' channel, and Ci and g are the LECs of the theory. Ellipsis 
denotes higher derivative four-meson terms which will be neglected in this section. Our goal is to determine the 
coupled-channel spectrum in a finite volume (FV) that falls within the p-regime of LQCD [45], [46], defined by 
3> 1 where ttItt is the pion mass, and L is the spatial extent of the volume. In this regime, finite volume 
corrections to the single particle dressed propagator are exponentially suppressed [19j . as are contributions from 
t,u-channel scattering diagrams [T^], [IHIIMI- The leading order (LO) volume effects in the two-body sector arise 
from the presence of poles in the s-channel scattering diagrams. These lead to power-law volume corrections to the 
2-particle spectrum |19j . |20j . Therefore in the following discussion we will restrict ourselves to the contribution of 
such diagrams to the scattering matrix A4, whose Mij matrix element corresponds to the scattering amplitude from 
the i*^ channel to the j*'' channel. Certainly, the t,u-channel diagrams contribute to the renormalization of the theory 



^ It is important to mention ttiat we are not making any claims about the relevance of the four-pion channel in the light-scalar sector 
of QCD. At this point, it is not evident how to incorporate such states into the calculation, therefore if one would choose to use the 
formalism presented here to study the light-scalar sector of QCD, there will be an overall systematic error that must bo accounted for 
through other means. 
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FIG. 1: Shown are the diagrammatic representations of the full scattering amplitudes. Solid (dashed) lines represent 
single particle propagators of the / (//) channel. The first two diagrams require the LECs to reproduce the 
scattering amplitude for the / channel. The grey circle (diamond) is the full scattering amplitude of the /(//) 
channel. The black dot denotes the effective coupling for the / channel, which includes an infinite series of 
intermediate // bubble diagrams. There are two other diagrams for the // channel that require the LECs to recover 
■Mil. II and have not shown. The third diagram ensures that the mixing term g reproduces the off-diagonal 
scattering matrix elements. 

and therefore to the definition of the LECs. But for energies bellow the inelastic threshold, their finite volume effects 
are suppressed. In the non-relativistic (NR) limit, the t,u-channel diagrams are prohibited, so at low energies it is 
sensible to neglect their contributions. Alternatively, one can redefine the LECs to absorb the contributions of the 
t,u-channel diagrams. 

When considering momentum-independent interactions, the infinite volume loops contributing to the s-channel 
diagrams are 
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2 J (27r)4 [{k - Py - m2 + ie] [k^ - + ie] ' 



(3) 



where P = {E, P) is the total four-momenta of the system. It is straightforward to write down the equations that 
{c/,c//,g} must simultaneously satisfy in order to reproduce the A4-matrix elements of the theory, similar to the 
approach of [17] in examining the role of AA intermediate states in the "^S'o NN scattering. These equations, which 
are diagrammatically shown in Fig. ([T]), can be written as follows 
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-M^,^ 



(4) 



where M.i^i is the full relativistic S-wave scattering amplitude for the i*'' channel, and A4jji is the amplitude 
responsible for mixing the two channels. Solving these coupled equations leads to a renormalized theory. 

Once this is done, one can study physics in a finite volume. In particular, we are interested in the energy-eigenvalues 
of the meson-meson system placed in a finite lattice with periodic boundary conditions. The spectrum can be 
determined by requiring the real part of the inverse of the finite volume scattering amplitude to vanish. In a periodic 
lattice, the integrals over the spatial momenta appearing in Feynman diagrams are replaced by a sum over discretized 
three-momenta. In finite volume the integral in Eq. ^ is replaced by 

dk" 1 

(5) 



riV — 



2L3 ^ 

k 



27r [{k - P)2 - m2 + ie] [k^ - m? + ie] ' 



where the spatial momenta are quantized due to the periodic boundary conditions k ~ 27rn/i for n e Z'^, while the 
temporal extent of our Minkowski space remains infinite. This sum suffers of the same UV divergence of Eq. (|3|, and 
the difference of the two, SGj = Gj — G°°, is finite. It is simplest to consider the case where the mixing term, g is 
small, and keep our expressions to leading order in g. Using this and the definitions of the LECs in Eq. Q, the finite 



5 



volume scattering amplitude of the / channel can be written as 

{Mij)v « 2 |,.2(.,+.„) ■ (6) 

Finally, we obtain the quantization condition for the coupled channel problem at LO in the mixing parameter 

ne[\Miji\'e^^(''+'"^ - (mij + (mhji + 1=0. (7) 



At this point we have refrained from using explicit expressions from the FV integrals and scattering matrix elements, 
these details will be presented in the following sections. It is important to note that despite the simplicity of this 
toy model, it illustrates all the features of the problem at hand, and, as it will be shown in the next section, the 
result for the weakly coupled channels presented above is the same to the strongly coupled case when only the s-wave 
contribution to the two-body scattering in the cubic lattice is taken into account. 



B. N arbitrarily strongly coupled two-body channels in a moving frame 

In the previous section, the following assumptions have been made. First of all, the higher order four-meson terms in 
the derivative expansion of the effective Lagrangian have been neglected. Secondly, the scattering is restricted to two 
coupled channels composed of identical particles, and only the s-wave scattering is considered. Most importantly, the 
mixing term was assumed to be small in deriving the quantization condition. As mentioned earlier, the latter condition 
is the most relevant when considering systems coupled via QCD interactions. In this section we will simultaneously 
remove all of this assumptions. 

Most of the details associated with generalizing to a moving frame have been developed by Kim et al. [22], [37) . 
which we will briefly review for completeness. A system with total energy and momentum, E and P, in the laboratory 
frame has a CM energy E* = ^/E'^ — P'^. For the z*'* channel with two mesons each having masses rrii i and mi_2, the 
CM relative momentum is 

which simplifies to ^ nif when rrii^i = 771^^2 = "^i- Including all possible momentum dependent operators in the 

EFT, the moving frame FV loop one needs to evaluate is [22], [S^ 




~ ^ y 27r [{k - Py - ml^ + If] [fc2 - rnf 2 + le] ~ ^ 

where Ui is 1/2 if the particles in the i*'' loop are identical and 1 otherwise. The function /(k) holds the momentum 
dependancc arising from the infinite tower of operators inserted on either side of the loop, i.e. the Bethe-Salpeter 
kernel, /C. The FV part of the sum can be written as 

SGY = - ^ £ E fLicuKM')] , (10) 

where /;,„ is defined via the spherical harmonic (i^m) decomposition of /(k*), 

CO / 

/(k*)=E E hramk*'^J^Y,^{^\<i>*). (11) 

The function c^^ySsfi^^ can be written in terms of the three-dimensional Zcta function, Zf^^^ 

C('^") = -^(t)'"^-^('^*^/2-)']' ZL{s:A-Y.^ (12) 

where the sum is performed over = {r e R'' | r = 7~^(m|| — ad) — m±,m £ Z}, d is the normalized boost vector 
d = PI//27r, the relativistic 7 is defined by 7 = E/E*, and a = ^ (l + with E* = \/ q* + m\ + q* +771^ 



dk° /(k) 
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The generalization of the quantization conditions for N channels that are coupled via an arbitrarily strong interaction 
is straightforward by upgrading the Bethe-Salpeter kernel, /C, to not just be a matrix over angular momentum but also 
over N channels. Similarly, it is convenient to represent the N loops as a diagonal matrix Q = diag{Qi,Q2, ■ ■ ■ ,Gn)- 
The kernel is assured to reproduce the infinite volume scattering matrix (M) by solving the following matrix equation 

iM - -ijc - ijcg°^jc - ijcg°°jcg°"jc + ■■■ = -nc- — I— ^jc = -m- — ^— -. (is) 

1 — y°°K, 1 — y°°J\/i 

With this definition of the kernel, one can proceed to evaluate poles of the N-channcls FV scattering matrix by 
replacing the infinite volume loops 0°° with their finite counterparts, 

^ ^{M)v = -^IC - ^ICg^lC ICG^ IC + ■■■ = -^IC^-^ = _ mG^)- (W) 

Finally arriving at the quantization condition 

Tie {det{l + 6G^ M)} ^ (15) 

where A4 is evaluated on-shell. For N=l, one reproduces the result first obtained by Rummukainen and Gottlieb [21] 
and later confirmed by Kim et al. [22] as follows. First note that it is convenient to evaluate the determinant using 
the spherical harmonic basis of SG^ and the on-shell scattering amplitude Aii [H] 

^ TL ' f 47r ^ 'sf f ^ 

{5GY)h,m^;h,m^ ^luh^m^m^ ~ * ~ XI 51 —^cf^ilf) / d^*Yil„,Yi*^Yi^„,^ (16) 

,,,, , , SttE* e^'^'P^it) - I 

{M^i)li,mx;l2,m2 = On 120 ml, m2 — . (17) 

If the two meson interpolating operator is in the irreducible representation of the cubic group, the energy 
eigenstates of the system have overlap with the I = 0,4,6, .. . angular states at zero total momentum, making the 
truncation at Imax = a rather reasonable approximation. When P ^ 0, the symmetry group is reduced, and at low 
energies the I = and / = 2 partial wave do mix Nevertheless, if one truncates the determinant over the angular 
momentum at Imax = 0, the familiar quantization condition 

4 cot(5°) = ^^^o'oli; i^L/^^n (18) 

is recovered. It is convenient to introduce a pseudo-phase defined by 

9*cot(0,) = 4^c^((Z*). (19) 

to rewrite the quantization condition as: 

cot((/)i) = — cot((5i) => 5i + (f)i — rmr, (20) 

where m is an integer. 

For the N=:2, the expression for the scattering amplitude in Eq. (17) is modified, as it now depends on the mixing 
angle e, and the scattering matrix is no longer diagonal, while still symmetric. The off-diagonal terms are labeled 
■Miji- Using the definition of the S-matrix for the coupled-channel system, Eq. Q, the scattering matrix elements 
can be written as 

(U ^ -A A 87r£;* cos(2e-)e^'^'"(^n-i 

UiQ* 2l 

{M^u)h,mv,i2,m2 = ^hM^mi.m^ , sin(2e) . (22) 

y^niniiq*jq}j 2 



Where the usual relativistic normalization of the states is used in evaluating the S-matrix element. From Eq. (15) 
one obtains 

T? ^ A + SGYMij SGYMiji \\ . 
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where the determinant is not only over the number of channels but also over angular momentum which is left implicit. 
In deriving this result we have made no assumption about the relative size between the scattering matrix elements, 
but when Imax = we recover the LO result in Eq. ([t]). For l„iax — one can show that the determinant has in fact 



no imaginary part, and using the pseudo-phase definition in Eq. (19), the quantization condition can be rewritten in 
a manifestly real form: 

cos 2e cos (01 + 5i - 02 - 52) = cos {4>i + (5i + 02 + <^2)- (24) 



One can see that in the e — ?► the decoupled quantization conditions for channels / and // , Eq. ( 20 ) are recovered. The 
extension to a larger number of coupled-channels is straightforward. As an example, we consider the N=3. Unitarity 
allows us to parametrize the S-matrix using three phases shifts {i5/, i5//, (5//} and three mixing angles {ei,e2,e3} 

/ e^'^^'ci ie'^^^'+^"hiC3 ie''^^'+^"''> S1S3 \ 

S3 = ie'(*^+*")siC2 e'25„ (g^g^gg _ ^2^^) i^riSr+S,ri) (g^g^Sg + S2C3) • (25) 

\ie^iSj+Sjir)s^s^ ie*('^"+'^"^) (C1S2C3 + C2S3) ie'^^'" (018283 - C2C3) J 

where q = cos(2ei), Si — sin(2ei). Note that in the limit £2 = £3 = the third channel decouples, and one obtains 
a block diagonal matrix composed of iS'2 corresponding to the I — II coupled channel, as well as a single element 
corresponding to the scattering in the uncoupled channel ///. The spectrum of three-coupled channel is defined by 

fl + SGYMij SGYMiji SGYMijii ' 

Tie { det SgYjMii.i 1 + SGYMuji SGYMujh \ } = 0, (26) 
\SGYhMihj dGYiiMinji I + SGYhMhi.iii , 



where the scattering matrix elements can be determined from Eq. ( 25 ) using the relationship between the scattering 
amplitudes and the S-matrix elements {Mi^j)i^^mi-M.m2 = &iiM^rni,m2^'^E* {{s''^^^)i^j - Sij) /{2iy^n^njq*q*). 



II. TWO-BODY ELECTROWEAK MATRIX ELEMENTS IN FV 



As discussed in the introduction, electroweak processes in the two-hadron sector of QCD encompass a variety of 
interesting physics, so it would be desirable to calculate the electroweak matrix elements directly from lattice QCD. 
One of the very first attempts to develop a formalism for such processes from a finite volume Euclidean calculation 
is due to Lellouch and Liischer. In their seminal work [33], they restricted themselves to fc — > tttt decay in the kaon's 
rest frame, and showed that the absolute value of the transition matrix element in an Euclidian FV is proportional 
to the physical transition matrix element. This proportionality factor is known as the LL-factor. This formalism was 
then generalized to moving frames by Kim et al. [55]. Here we present the generalization of Lellouch and Liischer 
formalism to processes where the initial and final states are composed of two-hadrons S-wave states. In the relativistic 
case, the coupled channels result Eq. ([t]) is used to derived the 2 — )■ 2 LL-factor for boosted systems. For the NR 
case, we use low energy EFT(7f) in two cases. The first case that is most pertinent for nuclear physics, regards the 
"^^o — ^Si mixing which requires including one-body contributions to the weak matrix. The relationship between the 
FV matrix element and the physical transition amplitude is then found. It has been shown that the absolute value of 
these two are no longer proportional to each other as a result of one-body weak contributions. We then also present 
the result for weak processes which do not involve any one-body contributions. 



A. Relativistic 2-Body LL-Factor 

1. Degenerate perturbation theory derivation 

In order to derive the relativistic two-body LL-factor, one first notes that in the absence of the weak interaction, 
the two states decouple, as a result the 5'-matrix becomes diagonal. As is pointed out by Lellouch and Liischer, there 
is a simple trick to obtain the desired relation by assuming the initial and final states to be nearly degenerate with 



energy Eq (each satisfying Eq. (20)) when there is no weak interactions. Once the perturbative weak interaction is 



turned on, the degeneracy is lifted, and the energy eigenvalues are 

E*^= E*±V\MYji\ = E*±AE, (27) 
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where n is the FV matrix element of the weak Hamiltonian density. Consequently, the CM momenta and the 
scattering phase shifts acquire perturbative corrections of the form 




V\MIjj\ ^ Aq* V\MYrj\ (28) 



P*3 
^0 



and 

ASM) ^ S'M)Aq* V\MYjj\. (29) 

where Sl{q*) denotes the derivative of the phase shift with respect to the momentum evaluated at the free CM 
momentum. The perturbed energy necessarily satisfies the quantization condition Eq. (23 1. The generalized LL- 
factor for 2—^2 scattering is then obtained by Taylor expanding Eq. ( 23 ) to leading order in the weak matrix element 
about the free energy solution, 

\MT;n\' = V \Aq*jAq*n (^) (^^) i^'M) + W^)) {<j>'Mqh) + S'jAqh))] IMjjA'- (30) 
I V niq*i J \niiq}jj J 

where (p'iiql) denotes the derivative of the pseudo-phase with respect to the momentum evaluated at the free CM 
momentum. 



2. Density of states derivation 



In the previous section, we arrived at the generalization of the LL factor for two-body matrix elements using the 
degeneracy of states argument. Lin et al. [51] showed that the the LL-factor for K ^ tttt can also be derived using 
the density of states in the large volume limit, and this argument was then generalized by Kim et al. |22j to boosted 



systems. Here it will be shown that the results in Eq. (30 1 is also consistent with the work by Kim et al. Let ai (x, t) 
be the two-particle annihilation operator for the z*'' channel. Then the two particle correlation function in FV can be 
written as 



C^, it) 



L— )-oo 



d^x e'P-^ (0| a, (x, t) a\ (0, 0) |0)^ = F ^ , 

m 

V J dEpvAE)e-''' m<j{O,0)\i;P,E)y\\ 



|cr(O,0) \i;P,m) 



(31) 



In the first equality we have inserted a complete set of states. In the second equality, we have introduced the density 
of states for the i*'' channel, pv,i{E), which is defined as pv.i{E) = dnii/dE. Using Eqs. (201, (28) the density of 
states can be written as pv,i{E*) ~ (0-(g*) + 5[{q*)) Aq* /tt. In infinite volume the two-particle correlation function 
is El] 



dE^e-'''ma{O,0)\^;P,E)^ 



(32) 



where the factor of has been introduced to account for the double counting of the phase space when the particles 
are identical. It is straightforward to show that this relation still holds when the two particles have different masses. 
From Eqs. (31 1, (32 1 the relationship between the states of infinite and asymptotically large (yet finite) volume can 
be deduced, 



l2Vpv,,E* 
n^q* 



V ■ 



(33) 



This relation therefore recovers the LL-factor as given in Eq. ( 30 1 . It also demonstrates that the LL-factor accounts 



for different normalizations of the states in the finite volume and infinite volume in presence of interaction. 



B. Two-Body LL-Factor in EFT(ji') 



In this section we discuss the generalization of the LL-factor for the two-nucleon sector. This sector has been 
previously studied by Detmold and Savage j40| in the finite volume. They considered a novel idea of studying 
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electroweak matrix elements using a background field. They point out that evaluating matrix elements of electroweak 
currents between NN states, e.g. < d\A'^\np >, is naively one or two orders of magnitude more difficult than 
performing NN-four point functions. This argument led them to present a procedure for extracting the relevant LECs 
of the pionless EFT, F,FT{ft) [15], by calculating four point functions of nucleons in a finite volume in the presence 
of a background electroweak field. This would be a project worth pursuing with great benefits, namely a five-point 
function is replaced by a four-point function, thereby dramatically reducing the number of propagator contractions. 
For isovector quantities, this procedure comes at a small cost, since for perturbatively small background fields the 
QCD generated gauge links get modified by a multiplicative factor that couples the valence quarks to the external field, 
U^'"^{x) — >■ U^'-^ ^ [x)U^^^ [x) . On the other hand, for isoscalar quantities this approach would require the generation 
of gauge configurations in the presence of the background field. For both isovector and isoscalar quantities, one would 
need to perform calculations at a range of background field strengths in order to precisely discern the contribution of 
the coupling between the background field and the baryonic currents to the NN spectrum. Additionally, the nature 
of this background field will defer depending on the physics one is interested in calculating. Alternatively, one can 
always evaluate matrix elements of electroweak currents with gauge configuration that solely depend on the QCD 
action, which is the case considered here. In the following section, we present the relationship between the physical 
transition amplitude for the isovector-axial current and its FV matrix element. As mentioned above, there are great 
benefits for performing this calculation using the background field approach, however it is necessary to have further 
crosschecks. 



1. ^Si) — ^Si mixing 

In this section, we will explore FV corrections of weak matrix elements in the two-nucleon sector. In particular, 
we will consider processes that mix the "^^o — '^5'i, e.g. proton-proton fusion — > d + e+ + i^e . In order to do this 
calculation the mechanism of EFT(;7f) [41], [42], [43] will be used. The methodology is similar to the one used for 
the toy problem considered in the first section, except the fields are now non-relativistic, and carry isospin and spin 
indices. The presence of a weak interactions, leads to a contribution to the Lagrangian that couples the axial-vector 
current A^^^ = ^ (uj^j^u — dj^j^d) to an external weak current. In terms of our low-energy degrees of freedom, the 
axial current will receive one-body and two-body contributions. At energies well bellow the pion-production threshold, 
the EFT(y) Lagrangian density including weak interactions can be written as [41], [42], [43] 

£ = iVt (^ido + ^^a\'^ N - Q (iV^P/iv) ^ (^N^P(n) (34) 

-C*o (N^PiNy (^N'^PiN^ - WsLi^A (iV^PfiV)^ (7V^P|7V) + h.c. +■■■ (35) 

where N is the nucleon annihilation operator with bare mass m, {Cq, Cp, ^i.a} are the LECs of the theory, 
gA — 1.26 is the nucleon axial charge, W3 is the external weak current, and {Pi, Pi} are the standard {^S'o,'^ 6*1}- 
projection operators. 

Pi = ;7|^2rV2 Pi = ^r^a.a^, (36) 



where T(cr) are the Pauli matrices which act in isospin(spin) space. In Eq. (34 1 the ellipsis denotes an infinite tower 
of higher order operators. The C'(p^")-operator for the {-^5o,^S'i} state will have a corresponding LEC {C'|„,C|„}, 
which are included in this calculation. In the absence of weak interactions, the scattering amplitude for both channels 
can be determined exactly in terms of their corresponding LECs by performing a geometric series over all the bubble 
diagrams [13] 

(2) Jj^w^^i^w:^, (3^) 

is linearly divergent. When defined with the power-divergence subtraction (PDS) scheme [HI [42l [52] , it becomes 

= -£ (a^ + ^VrnE - PV4) = £ (m + V) . (39) 



where the loop integral 



TOC 
^0 
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FIG. 2: Shown is the definition of the effective two-body couphng, ' (represented by the black circle), between 
the singlet and triplet di-baryon states. The grey circle denotes the bare two-body weak coupling and the crossed 
circle is the one-body weak coupling. 



where /i is the renormalization scale and as usual q* is the momentum of each nucleon in the CM frame. It is well 
known that the NR scattering amplitude can be parametrized in terms of the phase shift, S, and the CM momentum 

q* as 

^.^^l^^e-^^^S-l^ (40) 
mq* cot — i mq* 2i mq* 2i 

where we have used the fact that in the absence of weak interactions, the S'-matrix is diagonal and equal to e^"^^ . 
In the presence of the weak interaction, we use the "barred" parameterization in Eq. ([T]) to write the diagonal and 
off-diagonal scattering matrix elements 

4:TT Sii-1 47r cos2ee'^'^'so(3si) _ 1 

A^iSoC^Si) = — 77-- — ^ 77-- (41) 

^ ' mq* 2% mq* 2i 

47r Si2 47r sin2ee'(''%+'^'-'-=^i) , ds, +5,, ^ 

mq* 2i mq* 2 

In the absence of one-body contributions to the weak matrix elements, the determination of the QC for this system 
would be quite similar to our toy model in section |I A[ The mixing between the channels will lead to three coupled 
equations requiring the LECs to reproduce the infinite volume scattering matrix, thereby renormalizing the theory. 
In order to deal with the one-body operator, it is convenient to introduce a volume-dependent effective two-body 
coupling, defined in Fig. (|2|. The regularization condition for this system at LO in the weak interaction can 

be written as 

where we have introduced a volume dependent two-body coupling 

and the integral J^'^^^ is defined as J,^'*-^-* (E, P) = —OeIq^'^^^E, P). The axial charge is fixed at its physical value. 
Once the LECs are renormalized, one can move on to evaluate the scattering amplitude for each channel in the FV. 
The poles of the scattering amplitude satisfy 

I - M.SoMss.gAWSJ^f - (m.So + ^) [m^s, + j;^) } = (^5) 

where SJq = j}^ — and SIq = Iq — 1§° as before. In order to find the relation between 6Iq and the Liischer Zeta 
function, Zqq, it is convenient to add and subtract the infinite volume two-particle propagator at zero energy in Iq . 
Then evaluate one of them using PDS and the other using the a hard cut-off. A, leading to 

„N m 1 \ ^ 1 toA 

^ = -4.'^ + P ^ ¥^Er^ + 2^ 

k 47n m 



Ue 



therefore 



SJV^E, P) = ( ^Z^o [1; il*L/2nr] - M ^ M (-0) - ^) (47) 
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At leading order in the mixing angle e and the axial coupling gA Eq. ( 45 1 is in agreement with the result by Detmold 
et al. go]. 

With this, it is straightforward to obtain the relationship between the FV matrix elements of the Hamiltonian 
density and the physical matrix element. In the absence of weak interactions, the two dibaryon states are degenerate 
with energy Eq and momentum CM momentum Qq, satisfying the free quantization condition cot(— 0) = cot(^). As 
the weak interaction is turned on, the degeneracy is lifted, leading to a shift in energy equal to SE = V\M'^\, where 



|A^^| is the FV matrix element of the Hamiltonian density between the ^5*0 and '^S'l states. Expanding the Eq. (45) 
about the free energy, and keeping LO terms in the weak interaction, one obtains 

[\^^So-^s^-9AW,^J^^ - [^)\^' + Sk){<^' + S'.So)\MYs,-.sf- (48) 

This result shows that in order to determine weak matrix elements in the nucleon sector, it is not only necessary 
to determine the derivatives of the phases shifts in the ^5*0 and '^Si channels, but it is also necessary to determine 
the axial charge. There is no clear crosscheck for this result, since the density of states approach cannot reproduce 
this. The presence of the one-body operator makes the mixing between the two states non-trivial, therefore one would 
expect a more complicated relationship between the FV and infinite volume states than the one predicted via the 



density of states approach. It must be mentioned that Eq. ( 48 1 does not depend on the value of W3 since it factors 
out from both the FV and physical matrix elements. 



2. {BB)i -> {BB)i 



Lastly, we show how Eqs. 



(45), (48) can be extended to the low energy baryonic S-wave channels that couple 

NN and H 



exclusively via a weak tow-body interaction, which is pertinent for studying AA^ — > NN and H — > AA'' (with unphysical 
light-quark masses in order to assure the energy of the final states is bellow the pion production threshold) . For these 
energies, one can generalize the EFT(^) Lagrangian in Eq. (34) by replacing the nucleon field operator onto an SU(3) 
octet baryon field operator, B^, with flavor index i. We can loosen the restrictions on the mass, and let it depend on 
the baryon flavor. Besides the spin and isospin structure, this is quite similar to our toy model in section |I A[ For 
generality, the channels will be labeled by the subscripts {/,//}. It turns out that this general case can be obtained 
from Eq. (45) in the limit — 0, after performing the following replacements: 



2-d2 



2/ii = 



1 

vail 



1 

mi 2 



q* ^ q* = \ 2fi,E - 



(49) 



where rriij is the mass of the j particle in the r channel with CM momentum q* . With this, it is straight forward 
to show that the spectrum of this system satisfies 



UeUMi-iile 




(50) 



Eq. (24), 



Using Eq. (47), one can show that the term inside the braces above is in fact real, and can be rewritten exactly as 



cos 2e cos (0/ + (5/ - 0// - Sn) = cos (0// -t- Sjj + + 5ii). 



(51) 



The only difference between this result and Eq. (24) is the fact that for NR theories 7=1 and the dispersion 
relationship differs its the relativistic counterpart. 

Using the degenerate perturbation theory approach, the relationship between the FV matrix elements of the Hamil- 
tonian density and the physical matrix element is found to be 



\MT_ 



(52) 



In this case, this NR result for the LL-factor can be confirmed using the density of states approach. The FV four-point 
correlation function, Eq. (31 1, is the same for relativistic and NR theories. On the other hand, in the infinite volume 
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the phase space and the dispersion relation differ between the relativistic and non-relativistic theories. That being 
said, the NR four-point correlation function in infinite volume is 



C 



|a(O,0)|z;P,ii;) 



NR,oo 



(53) 



Comparing this expression of the infinite volume correlation function to that of the asymptotically large FV results, 



Eq. (31), we find the relationship between the FV and infinite volume NR states 



\i;F,E) 



NR.oo 



4^ 



V,i 



NR,V ' 



PV,i 



which is consistent with Eq. ( 52 ) 



471(7* ' 



(54) 



III. SUMMARY AND CONCLUSION 



In this paper, we have presented and derived the FV quantization condition for a system of multi-coupled channels 
each composed of two-hadrons in a moving frame. In the first section, the quantization condition at LO in mixing 
phase, Eq.([7]) is derived for the s-wave scattering, where a rather simple EFT toy model for scalar particles is used. 
Using the techniques developed by Kim et al. |22) . the quantization condition for boosted systems with arbitrary 



mixing angles, Eq.( 15 ) has been obtained. It then became evident that the toy model used is in perfect agreement with 
the non-perturbative result when the series over angular momenta is truncated at the s-wave. From the generalized 
results, we have also derived the quantization condition for N=3 coupled channels composed of two hadrons, Eq.(26|. 



This result has the advantage that it would allow the lattice practitioners to perform coupled channel calculations 
at multiple boosted momenta in a periodic lattice, thereby increasing the number of measurements in order to best 
constrain the S-matrix elements. When N=2, the S-matrix can be parametrized by three real parameters (two 
scattering phase shifts and one mixing angle) , therefore one needs to perform at least three measurements at each CM 
momentum, which can be done by using combinations of different boost momenta and different volumes. The N=3 
case would require six measurements for each CM momentum to constrain the three phase shifts and three mixing 
angles. 

We have also derived the relationship between FV matrix elements and infinite volume matrix elements in the 
two-body sector, Eq. ( 30 ) . This was first calculated for K ^ tttt in the rest frame by Lellouch and Liischer [39] and 
later extended to the boosted system by Kim et al. Here we showed two ways to obtain the extension of LL-factor 
for 2 — )• 2 relativistic processes. The first entails expanding the coupled quantization, Eq.Q, about the free energy 
of the system when the two channels are decoupled. This method assumes the two states to be degenerate in the 
absence of the weak interaction, as was first developed by Lellouch and Liischer [39]. Kim et al. [22] generalized the 
method of |51j and derived the relationship between FV and infinite volume two particle states in the moving frame. 



Eq. ( 33 1 , which has been shown to agree with our result of the two-body relativistic LL-factor Eq. ( 30 1 



We have used EFT('ji') [1TH151 [551 to determine the extension of the LL-factor for non-relativistic baryonic 
systems. First processes that mix the ^Sq — ^Si NN channels have been considered. This is pertinent for performing 
calculations of proton-proton fusion, among other interesting processes, directly from LQCD, [ID]. The channels in this 
system are mixed not only by a two-body operator but also by a one-body operator. As it is shown, FV effects arising 
from the insertion of a one-body operator are sizable and therefore must be included. It is convenient to introduce a 
volume-dependent effective mixing vertex that encompasses contribution from the one-body and two-body operators, 
Eq.(44|. The effective vertex is required to reproduce the infinite volume transition amplitude. Expanding our 



quantization condition, Eq.(45l, about the free energy solution, we find the relationship between the FV and infinite 
volume weak matrix elements, Eq.(48|. Unlike any previous case, the FV and infinite volume weak matrix elements 
are not simply proportional to each other, and their relationship depends on the nucleon axial charge. We do not 
present a crosscheck for Eq.(48), because it is not clear how to extend the density of states approach in the presence 
of one-body operators. Eq.(48) allows one to determine physical weak matrix elements in the two-nucleon sector from 



LQCD, up to exponentially suppressed finite volume effects. 

Finally, we have derived the NR coupled channel quantization condition for systems that do not have a one-body 



current operator contributing to the physics of interest, Eq.(50), which would be relevant for studying the H-dibaryon 
decays. This result is universal, and is independent of the spin/flavor composition of the two particle states. As in 
the previous case, the NR coupled-channels quantization condition can be expanded about the free energy solution to 
derive the LL-factor for 2 — >■ 2 NR processes in the absence of one-body contributions to the weak matrix elements. 



13 



Eq.(52|. This result allows the practitioner to determine the dominant FV contributions arising when evaluating 
matrix elements of electroweak currents. In order to confirm this result, we derived the relation between the FV and 
infinite volume two particle state for NR systems using the density of states approach, Eq. ( 54 1 . We found agreement 
between the two procedures, further confirming Eq.(52|. 
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